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Numerical investigation of the quantum fluctuations of optical fields transmitted 

through an atomic medium 
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We have numerically solved the Heisenberg-Langevin equations describing the propagation of 
quantized fields through an optically thick sample of atoms. Two orthogonal polarization compo- 
nents are considered for the field and the complete Zeeman sublevel structure of the atomic transition 
is taken into account. Quantum fluctuations of atomic operators are included through appropri- 
ate Langevin forces. We have considered an incident field in a linearly polarized coherent state 
(driving field) and vacuum in the perpendicular polarization and calculated the noise spectra of the 
amplitude and phase quadratures of the output field for two orthogonal polarizations. We analyze 
different configurations depending on the total angular momentum of the ground and excited atomic 
states. We examine the generation of squeezing for the driving field polarization component and 
vacuum squeezing of the orthogonal polarization. Entanglement of orthogonally polarized modes is 
predicted. Noise spectral features specific of (Zeeman) multi-level configurations are identified. 

PACS numbers: 42.50.Ct 42.50.Dv 42.50.Gy 03.67.Mn 42.50.Nn 
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I. INTRODUCTION 

The preparation of optical fields in states with purely 
quantum-mechanical properties is the key ingredient of 
quantum optics and the essential requirement for their 
use in quantum information processing. Light fields 
presenting squeezing are well known examples of non- 
classical states for which numerous applications have 
been suggested and demonstrated. Entanglement of a 
two-mode field is another important example of a purely 
quantum mechanical resource which lies at the basis of a 
large number of quantum information procedures such as 
EPR pair production, telcportation, and quantum cryp- 
tography For two degenerate field modes, considered 
as continuous variable systems, squeezing and entangle- 
ment are related concepts. Squeezing in one mode leads 
to entanglement between two modes obtained by a lin- 
ear optical transformation, as implemented by a beam 
splitter HQ. 

Ever since the first proposals for light squeezing, 
atomic systems have attracted considerable attention ow- 
ing to the large nonlinearities associated to resonant 
transitions. Indeed, the first successful demonstration 
of squeezing used non-degenerate four-wave mixing in 
sodium atoms [J] contained in an optical cavity. Several 
subsequent experiments 0, 0] also used atomic samples 
contained in optical cavities. In these experiments, the 
highly nonlinear interaction between the atoms and the 
cavity-mode plays an essential role in the generation of 
squeezing. 

Following the early work, most present day exper- 
iments on the generation of non-classical fields with 
atomic samples involve the use of optical cavities and 
require rather complicated experimental setups [7], la, |9| • 
However, in view of applications, the use of single path 
schemes for the generation of non-classical light fields 



could be of considerable practical interest. In this pa- 
per, we are concerned with the modification of quantum 
fluctuations of a single monochromatic light beam inter- 
acting with an atomic medium on a single path. Such 
possibility is alre ady implicit in the pioneering work by 
Walls and Zoller [l(| and Mandel [III predicting reduced 
quantum fluctuations in the light emitted by a resonantly 
driven two-level atom. The spectral distribution of the 
quadrature fluctuations of light emitted by a driven two- 
level atom was first calculated in [l2j]. The generaliza- 
tion of this study to an extended atomic sample was 
carried out by Heidmann et al. [Hj], who considered 
the fluctuations in the field emitted in the forward di- 
rection by a thin layer of atoms at rest, driven by a 
normally incident laser beam. At low laser intensities, 
squeezing is predicted for the low frequency components 
of the in-phase quadrature. For saturating intensities, 
squeezing occurs for noise frequencies around the gen- 
eralized Rabi frequency £1 = (A 2 + SI 2 ,) 1 / 2 , where A is 
the atom-laser detuning, and flo the incident field reso- 
nant Rabi frequency. The calculation of the fluctuation 
spectra of light traversing a thick two-level atom medium 
was presented by Ho et al. Il4i. Single path squeez- 
ing was observed using sodium[15| and ytterbium 
Quite recently, intensity-intensity quantum correlations 
(sub-shot noise intensity-difference noise) were observed 
in non-degenerate forward four-wave mixing in rubidium 
vapor [l7| . 

The initial work on squeezing through light-atom inter- 
action considered ideal two-level transitions and a single- 
mode optical field. A more realistic approach requires 
the consideration of multi-mode light fields (including 
different polarizations) and multi-level atoms. Three- 
level atoms interacting with two fields (A system) have 
been analyzed. Such systems bear a considerable inter- 
est owing to the possibility of large nonlinearities in as- 
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sociation with electromagnetically induced transparency 
(EIT). Phase- noise squeezing was predicted for a A sys- 
tem in which one of the fields was classical, for a cavity 
contained atomic system (l8j and for single path propa- 
gation through an ensemble of motionless atoms [19( . 

An interesting issue of the multi-mode field interaction 
with an atomic sample is the possibility to achieve po- 
larization squeezing, which is signalled by squeezing of 
the vacuum field with orthogonal polarization relative to 
the incident field. Polarization squeezing with cold atoms 
inside an optical cavity was demonstrated by Josse and 
co-workers [1, |2(J. Polarization squeezing is intimately 
related to the observation of continuous variable entan- 
glement between two field modes [H, [2l|, [H, HH . 

Recently, a renewal of attention on the fluctuations 
of light transmitted through an atomic sample was moti- 
vated by the prediction that vacuum field squeezing could 
be achieved on a single path as a consequence of polar- 
ization self-rotation (PSR) [24|. An experimental obser- 
vation of squeezing via PSR was reported [25[, using a 
room temperature rubidium cell. Such result could not 
be reproduced by other groups, in spite of the method- 
ical exploration of the relevant experimental parameter 
space 26]. It is argued by Hsu and co-workers (2(| that 
excess noise, preventing the observation of squeezing via 
PSR, originates from the quantum atomic fluctuations 
not explicitly included in the original theoretical pro- 
posal [24j . Such argument is supported by a simplified 
four-level system model 20], for which the noise arising 
from atomic quantum fluctuations dominates over the 
semiclassical squeezing terms, under the conditions cor- 
responding to the experiments. However, PSR squeezing 
is not excluded for cold atom samples in the regime of 
large intensities and optical detunings. 

The noise properties of light fields interacting with 
atom samples on a single path have been experimen- 
tally investigated under conditions of EIT. Large corre- 
lations and anti-correlations were observed between the 
two fields at the Raman resonance condition between the 
two ground state hyperfine levels of Rb [13, [28[. The 
change in sign of the correlation is determined by the 
light intensity. Noise spectra and correlations between 
two polarization components participating in Hanlc/EIT 
resonance involving Zeeman sublevels of the Rb ground 
state were observed [29| . These studies were carried out 
using diode lasers, known to possess large excess phase 
noise. A qualitative agreement between these experi- 
ments and theory could be reached by considering a clas- 
sical field with random phase diffusion. 

All theoretical models considered so far in connection 
with the analysis of light fluctuations interacting with 
an atomic sample are based on several simplifications. 
Single-mode approaches ignore the light polarization or- 
thogonal to the incident field, through which vacuum 
fluctuations enter the atomic system. Such fluctuations 
interact with the incident field, provided the sample pos- 
sesses or acquires some anisotropy through interaction 
with light. Multi-mode models also rely on simplify- 



ing assumptions. Sometimes one of the incident fields 
is taken as classical [llj]. The transverse spatial struc- 
ture of the field is generally not considered [3(| • In most 
cases, quantum fluctuations of the atomic operators are 
ignored. Even when a full quantum treatment was used, 
the atomic level structure considered was assumed to be 
an ideal three or four-level system. In most studies, the 
effect of propagation through an optically thick sample is 
not examined. The effect of the atomic velocity distribu- 
tion and that of nearby atomic levels within the hyperfine 
structure is usually neglected. 

It is the purpose of this work to investigate the quan- 
tum noise in fields transmitted through an homogeneous 
atomic sample when the full Zeeman degeneracy of the 
ground and excited atomic levels is taken into account. 
We consider a light field incident on the atomic sample 
with a well defined polarization and take into account 
the incoming vacuum field fluctuations with orthogonal 
polarization. We calculate the field fluctuation noise 
spectrum after propagation through the medium for 
arbitrary polarization and field quadrature angle, fully 
taking into account the influence of the quantum 
fluctuations of the atomic medium. Depending on the 
choice of the ground and excited level angular momenta 
(F g and F e respectively) and the field polarization, 
several configurations can be analyzed. A two-level 
system is obtained for a F g = — > F e = 1 transition. 
An open A system is obtained for a F g = 1 — ► F e = 
transition, if the incident field is linearly polarized and 
the two circular polarization components are considered. 
The four level system studied in [2(| corresponds to 
a. F g = 1/2 — > F e = 1/2 transition, with linear field 
polarization. In addition to those schemes previously 
explored, our calculation allows us to address other 
configurations, such as the F g ^ — > F e = F g + 1, 
for which no full quantum treatment was previously 
reported in spite of the prediction of PSR squeezing 
in this system [24[. Our calculation considers a 1- 
dimensional propagation (along the z axis) through a 
spatially homogeneous atomic sample considered, as a 
continuous medium [3l| . We begin by assuming that the 
atoms are at rest. Next, we briefly address the effect 
of a thermal atomic velocity distribution. Propagation 
effects on quantum fluctuations for a thick medium, of 
length L, are included in the calculation, although the 
depletion of the incident field mean value is ignored for 
simplicity. Thus, our results can be directly applied to 
situations in which absorption of the carrier frequency 
component of the field can be neglected (saturating field 
intensities or large atom-field detunings) . 

Our calculation confirms the prediction of quadra- 
ture squeezing under several conditions and in- 
dicates that squeezing by PSR can take place 
for F 9 = 1/2 — > F e = 1/2 transitions, and for 
F g > — > F e = F g + 1 transitions as well. The rather 
complex structure of the noise spectrum for thick optical 
media is illustrated and new features in the noise 
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FIG. 1: (Color online) Scheme of the physical situation the- 
oretically addressed in this paper (Pi, P2: polarizers, W\, 
W2: waveplates). The driving field polarization is imposed 
by Pi and W\, while W2 and P2 achieve the decomposition 
of the transmitted field into orthogonal polarization compo- 
nents. The noise analysis is made by homodyne detection. In 
the calculations, we have chosen linear x polarization for the 
driving field and analyzed fluctuations for polarizations x and 

y- 



spectrum specific of the multi-level Zeeman structure 
are presented. 

The paper is organized as follows. In the next section, 
the main lines of our theoretical model and noise 
spectrum calculation are presented. In section IIII1 we 
present and discuss the noise spectra corresponding to 
different level schemes. Concluding remarks are made in 
section HVl 



II. MODEL 

The problem that we address in this paper is schemat- 
ically presented in Fig. [TJ A laser beam traverses an 
optically thick atomic sample. The polarization of the 
incident driving field is defined by appropriate polariza- 
tion optics. After the sample, the field is decomposed 
into two chosen orthogonal polarizations. The quadra- 
ture noise on both polarization components is analyzed 
by homodyne detection [32|. 

Our theoretical approach follows closely the method 
used by Dantan and co-workers [33j for the study of the 
propagation of light fluctuations through an ensemble of 
three-level atoms in a A configuration. We consider two 
levels with Zeeman sublevels: a ground state g, of total 
angular momentum F g and zero energy, and an excited 
state e, of angular momentum F e and energy Hujq. The 
total radiative relaxation coefficient of level e is T. We 
assume that the atoms in the excited state can radia- 
tively decay into the ground state g at a rate bT, where b 
is a branching ratio coefficient that depends on the spe- 
cific atomic transition (0 < b < 1). For a closed (cycling) 
transition, 6=1. When the transition is open (b < 1), 
excited atoms can decay back into level g or into levels 
external to the two-level system. The atoms are in the 
presence of a magnetic field B directed along the light 
propagation axis z. In order to simulate the effect of a 
finite interaction time of the atoms with the light, we 



introduce an overall phenomenological decay constant 7, 
(7 <C r) which is compensated, in the steady-state, by 
the arrival of fresh atoms in the ground state. We ini- 
tially analyze an homogeneous ensemble of atoms at rest, 
leaving the effect of the atomic velocity distribution for 
subsequent consideration. 

Incident upon the atoms is a light field described in 
the Heisenberg picture by the operator: 

E(z,t) = £,(a 1 e i{kz -" Lt) e* 1 +a 2 e l{kz - u>Lt) e* 2 (1) 
+a\e- l{kz - u>Lt) e l + a^ 1 ^^^) , 



where e\ 
larization 



and 62 are two orthogonal (complex) po- 



unit vectors and 01, a 2 , 



the 



slowly varying field operators obeying the commutation 
rules [a K (z, i), a\[z' , t')\ = and a K (z,t),a\(z' ,t') — 

5 K \^5{t — t' — The quantization length L is cho- 

sen as the atomic medium length and c is the speed of 

light in vacuum. £ = <J 2 e "si, 1S ^ ne sm 8'l c photon field 
amplitude (where S is the mode cross-section). 

The atomic operators in the Heisenberg representa- 
tion for an atom j at position Zj are p^„(i) = 
where |^) and \v) designate Zeeman sub-states. We in- 
troduce the slowly varying atomic operators cr^(zj, t) = 
U^ v {t)W where U = p3 e ^-^ L t) + p j is an unitary 

transformation. P-j and P| are the projectors on ground 
and excited state manifolds respectively. Following [33| . 
we define continuous local operators (at position z) by 
averaging over a slice of the atomic medium of length 
Az: 



a^(z,t) = lim 



Az^o NAz 



E 



( Zj ,t), (2) 



z^Zj ^z+Az 



where L is the total length of the atomic medium and N 
the number of atoms. 

The atomic Hamiltonian is: 



n r 

H A = — (H +H B )dz , 



(3) 



with Hq = hwfjP e the isolated atom Hamiltonian and 
Hb = {PgPg + PePe)FzB the Zeeman coupling with the 
magnetic field B. P e and P g are the local projectors on 
the excited and ground manifolds, respectively, (3 g and (3 e 
are the ground and excited state gyromagnetic factors, 
and F z is the local total angular momentum operator 
projection along the magnetic field axis z. 

The atom field coupling Hj ntl in the rotating wave 
approximation, is: 



H 



Int 



~L 



(a\ei + c4e2). 



h.c. 



dz 



(4) 



Here Q ge — {QegY — PgQPe is a dimensionless op- 
erator related to the atomic electric dipole operator D 
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through D = (g\\D\\e) Q. The reduced matrix element 



<?||.D||e ) of the dipole operator between the ground and 



excited state is taken to be real, i] = £ (g\\D\\e) jh is 

the reduced atom-field coupling constant (half the single 
photon Rabi frequency). The standard spherical compo- 
nents of the operator Q eg are Q\ gl with (q = —1,0, 1). 
Their matrix elements Q q ge = \ Q q ge \v) (where \x 
and v refer to Zeeman substates belonging to the ground 
and excited manifolds respectively) are the correspond- 
ing Clebsch-Gordan coefficients. 

The complete set of atomic operators can be organized 
into a two-dimension operator array a = {cry} whose el- 
ements are the individual oij operators. After manip- 
ulation, one can formally write the Heisenberg-Langevin 
equations (including relaxation terms) for the atomic op- 
erators in the form: 



The mean value of the atomic operators, given the in- 
cident field mean values, can be obtained by taking the 
mean value of both sides of Eq. O This gives the usual 
Bloch equations. A numerical solution of these equations 
was presented in [34j. We consider next the evolution of 
the fluctuation operators up to first order. Fourier trans- 
forming the first order contributions in Eqs. [5] and [5] we 
get: 



d_ 

dz 



—S aj (z,w) 



— Sa\(z,w) 
oz J 



i—5aj(z, uo) 



(7a) 



+iN-e.j. ^ Qge,ikSa k i(z,uj) 



= i—5aUz,uj) 



(7b) 



iN^e*. S<j lk {z,uj)Q 



eg,kl 



da/dt 



-iA [P e ,a]--[H B ,a] 



+irj 



|(aie* 



a 2 e* 2 ).Q e 



h.c 



-bT(2F e + l)} ^QlaQl g 



9 



.{P e ,cr}- 7 (<7-Cr )+/ 



(5) 



where A = ljq — lul is the optical detuning and 700 
is a pumping term describing the isotropic arrival of 
fresh atoms in the lower ground state, er° = {cr^ } is 
chosen such that <r°- = for i ^ j, 0% = ( 2 f I +i) 
(I is the identity operator) if state i belongs to the 
ground level, and er° = otherwise. / = {/y} repre- 
sents the set of Langevin force operators satisfying [33| : 

(fv(z,t)fti( z ',t')) = j?2D ijtia 6(z'-z")8(t-t) , where 
Dij t ki is the corresponding diffusion coefficient. 

The field evolution is governed by the Maxwell- 
Heisenberg equations (in the slowly varying envelope ap- 
proximation) : 



.d d . 
{ dt +C d-z )aX 



,d d f 
■di +C d-z )a " 



,ik&ki (6a) 

-iNl]i* x . ^ VlkQegM (6b) 



with A = 1,2. 

Eqs. 03 and [5] form a set of coupled differential equa- 
tions for the atom-field interaction along the atomic sam- 
ple. We simplify the solution of these equations by 
assuming that the process is stationary and that the 
incident field mean value is undepleted as the atomic 
medium is traversed (a\(z)) = (a\(0)) = (a\). We 
proceed by linearizing the field and atom operators: 
(Tij(z,t) = (aij)+5aij(z,t) and a\(z,t) = (a x ) + 6a\(z,t) 
with (5(Jij) = (Sa\) = 0. 



jS<7 



-IT] 

-ii] 



iA[P e ,5a] + -[H B ,5a] 

h 



|((<4) ei + ^4) C2)-Qge + h.C.} ,5(7 

^(Sa\ii + 5a\e 2 ).Q ge + h.c.} , (a) 



-bT(2F e + l)(J2Q q ge ScrQl g ) 



--{P e ,5a} + "/8a 



(8) 



where f(uj) 
^2D ijtk iS(z 



i {f i3 (z,uj)}, with ^i(2,w)/i(z'w')| = 
z')5{uj-u)'). 

Field fluctuations depend linearly on atomic fluctua- 
tions, which in turn are driven by the Langevin force op- 
erators. To numerically solve these equations, we adopt 
a Liouville-space approach [34l [35j , organizing all opera- 
tors Uij into a column vector x, with n — 4(F g + F e + l) 2 

elements, and the four field operators ai,a\,a,2,a 2 into 
a four-element column vector A. Then, with some ma- 
nipulation, Eqs. [7] and [8] can be written in the form (we 
drop the dependence on z and lo for shortness): 



dSA .oj v . . Ntj 

-— = t-l 4 SA H -WSx 

oz c c 



- {iuA n + A)8x = f + r/V6A 



(9) 



(10) 



where I„ is a n x n identity matrix, W is a 4 x n matrix 
dependent on the coefficients Q?- , A is a n x n matrix cor- 
responding to the atomic evolution (including relaxation 
terms), and V is a n x 4 matrix describing the coupling of 
the field fluctuations to the atomic operator mean value. 
By defining M = — (iojI n + A), we can invert Eq. [TO] and 
from Eq. [9] we get: 



oz c 



(11) 
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withG = WM- 1 a 4xn matrix and B = {i^h+^-GV) 
a 4 x 4 matrix. A formal solution of Eq. [Tl] for propaga- 
tion over a length z is given by: 



6A(z,w) 



e Bz [6A(0,u) 



(12) 



Nr, 



-Bz' 



Gf{z',cu)dz' 



The power spectra of the fields' fluctuations after prop- 
agation through an atomic medium of thickness z can be 
obtained from the matrix S(z,u>) related to the field op- 
erators' spectral correlation matrix through: 



light fluctuations owing to the mean value of the 
atomic polarization in response to the mean incident 
field. Such a term may lead to noise reduction and 
cross polarization effects. The second term on the 
right-hand side of Eq. [18] represents the light noise 
introduced by the atomic quantum fluctuations via the 
matrix X, which is determined by the Langevin forces' 
diffusion matrix D (Eq. I16p . It always corresponds 
to noise increase (detrimental to squeezing). The 
calculation of the diffusion matrix D can be made with 
the help of the generalized Einstein theorem [37| and 
the corresponding expressions are given in the Appendix. 



(SA(z, u) [5A(z, io')} f ) = ^S(z, u)S(lu - u 1 ) . (13) 

Making use of (f(z, uj)P{z', a/)) = ^2D6(z-z r )S(uj- 
a/), where D is the atomic Langevin forces' diffusion ma- 
trix, and of Eq. [T2]one obtains: 



S(L,cj) = e BL S(0)(e BL y + 2^e BL x 



>~ Bz ' GDG^ e~ Bi z ' dz' 



(14) 



where we have omitted the dependence on ui for shortness 
and taken z = L. 

The term proportional to the identity in B commutes 
with all other operators, so we can write: 



S(L,u) 



KL afn\„K^ L , Nr, K L 



s(oy 

L 



e- K *'je 



e 

-A't 2 



(15) 



dz' 



In order to be able to calculate noise spectra for ar- 
bitrary field quadratures, we introduce the rotated field 



operators array 8A(9) = | 



o,\e , a^e , die , a 2 e 



'}■ 



It can be immediately calculated as SA(9) = T 5A, 



where T is the 4x4 matrix with e 



-w 



-i0 J.0 



along 



its main diagonal. In a similar way, fluctuations along 
two arbitrary orthogonal polarization unit vectors e[ and 
e! 2 can be calculated as SA'(9) = R SA(6), where R is 
the polarization-basis-change matrix. Thus, the spectral 
density matrix S' g {L, ui) for a given quadrature angle and 
arbitrary choice of the polarization basis can be evaluated 
as: 



S' 9 {L,lu) =TRS{L 1 Lu)R jf T j < 



(19) 



The noise spectrum s'-g for the quadrature angle 9 of the 
field with polarization e'j (j = 1,2), as can be measured 
by homodyne detection, can be evaluated from the ma- 
trix elements S' e of S' g (L,ui). Dropping, for shortness, 
the L and u) dependence we have: 



where J = 2GL>G t and K = ^-GV. 

c 

Let X' be a matrix satisfying: 

- (KX' + X'K ] ) = J 



(16) 



Then the integral in Eq. [15] can be evaluated and we get: 



S(L) = e KL S{0)e K L 
Nr, 2 



(17) 



(X' - e KL X'e K " L ) 



Introducing X — E^- we have: 



S(L) 



e crGV S(0)e 



(ctgv) 1 ' 



(18) 



-C 



X - e CTGV Xe {CTGV)l 



where C — N, ( ! r L is the cooperativity parameter [33L |36] | . 

Eq. [TBI gives the spectral density matrix S(L,uj) after 
the atomic medium, given the incident field spectral 
density S(0,lj). The first term on the right-hand 
side of Eq. [TS] represents the semiclassical effect on 



022 



^2(1 



934 



^043 



S'n 



044 



(20a) 
(20b) 



Other matrix elements of S g not appearing in Eqs. [20] 
are related to the field correlations between the two po- 
larizations. 



A. Velocity distribution 

So far we have considered an homogeneous sample of 
atoms at rest. We will generalize now our calculation to a 
sample of moving atoms with velocity v z in the direction 
of the light propagation axis. The velocity distribution 



is W (v z ) obeying /_ 
becomes: 



+ OC 



W(v z )dv z = 1. Then Eq. [TT] 



d[5A(z,uj)} _ B , 5A | Nr, 
dz c 



W(v z )G(v z )f v Jv z , 
(21) 
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with 



B' = 



i-I 4 

c 

Nrf 



(22) 



+ 00 



W(v z )G(v z )V(v z )dv z 



G(v z ) 
M(v z ) 



WM(v z ) 1 
—twin + A(v z ) 



(23) 
(24) 



A(v z ) depends on v z through the velocity dependent 
detuning A = u>q — lol — kv z (see Eq. [5]) and V (v z ) 
depends on v z through the mean values (x (v z )). 



After formal integration of Eq. [21] one has: 



5A(z,u>) 



„B'z 



Nrj 



-B'z' 



(25) 



+ OC 



W(v z )G(v z )f(v z ,z',uj)dv z dz' 



Since the Langevin forces for atoms with different ve- 
locities are uncorrelated, we impose: 



(f(z,v 2 ,oj) /t (*»')) 



(26) 



L 



NW (v z ) 



2D5 (z - z') 5{uo- uj') 6 (v z - v' z ) 



Using Eqs. l25l and l26l one can obtain the spectral den- 
sity matrix for the output field as: 



S{L) 



with 



- K ' z 'j'e- K " z 'dz' 



(27) 



and 



K'L = CT 



J' = 



+ 00 



W(v z )G(v z )V(v z )dv z , (28) 



+ 00 



W{v' z )G{v' z )DG{v' z )Uv[ 



(29) 



The expression for the spectral density matrix given in 
Eq. [27] represents a generalization of Eq. [15] that can be 
evaluated by the same method. The integrals in Eqs. [25] 
and [2H can be evaluated numerically, given the velocity 
distribution W(v z ). 



III. CALCULATED NOISE SPECTRA 

The model developed in Section [IT] allows us to calcu- 
late the spectral density of the fields transmitted through 
an homogeneous atomic sample of thickness L and total 



number of atoms N under a wide range of conditions, 
i.e. arbitrary choices of the atomic levels' angular mo- 
menta, driving field polarization, and incident field mean 
intensity. In addition, the model allows the arbitrary 
choice of the branching ratio b and the longitudinal mag- 
netic field B. However, in this paper we will restrict 
ourselves to consider a reduced number of parameters. 
No magnetic field is considered and the branching ratio 
b is taken as one (no transition to states external to the 
two- level system). Nevertheless, one should notice that 
some choices of F g and F e effectively describe open sys- 
tems since there are trapping (dark) states not coupled 
to the applied field. The driving field polarization (ei) 
is chosen to be linear along axis x, which means that 
vacuum enters the system with y polarization (£2). The 
x polarized driving field is assumed to be in a coherent 
state I a). The corresponding reduced Rabi frequency (in- 
dependent of polarization and atomic dipole orientation) 
is il r = 2arj, taken real. Under these conditions, the 
incident (white noise) spectral density matrix is [351 ] : 



S(0,w) 



10 



10 





(30) 



in units for which the standard quantum noise limit cor- 
responds to one. 

We have analyzed the output field fluctuations for the 
same set of orthogonal polarizations [x and y) used for 
the incident field (R = I4 in Eq. [H?]) . Noise spectra for 
a different choice of the polarization basis can be easily 
obtained using the appropriate matrix R. We have cal- 
culated the spectra for both the amplitude quadrature 
(6 = 0) and the phase quadrature (6 = tt/2). In all cal- 
culations, the relaxation rate 7 was taken as 7 = 0.0 IT, 
a realistic figure for experiments using alkali atoms' D 
lines. In the following, we analyze different level schemes 
depending on the angular momentum of the ground and 
excited states. We begin by considering the case of an 
homogenous sample of atoms at rest (cold atom sample) . 
The effect of atomic motion is discussed at the end of 
this section. 



A. F g = -> F e = 1 

Choosing the quantization axis in the direction (x) 
of the linear polarization of the applied field, this case 
reduces to a pure two-level system in which the field 



couples the Zeeman sublevels \F„ 



0, m a 



0) and 



\F e = l,m e = 0). The y polarization, along which vac- 
uum is entering the atomic system, is totally uncoupled 
to the driven transition. As a consequence, no change 
takes place for the y polarization mode and this field 
emerges from the sample in the vacuum state. 

The results of our model are in agreement with pre- 
vious work on quantum fluctuations and squeezing for 
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FIG. 2: (Color online) Transition F a = -> F e = 1 (Inset: 
the level scheme is shown using x as quantization axis; hollow 
arrow: driving field). Noise spectra of the amplitude (solid) 
and phase (dashed) quadratures of the transmitted field with 
the same linear polarization as the driving field, at zero de- 
tuning. The different values of the reduced Rabi frequency 
fl r (in units of T) are indicated (C = 1, 7 = O.Oir). 



FIG. 3: (Color online) Transition F g = -> F e = 1. Noise 
spectra of the amplitude (solid) and phase (dashed) quadra- 
tures of the field transmitted with the same linear polariza- 
tion as the driving field for A = 10r = fi r = 10r. The dif- 
ferent values of the cooperativity parameter C are indicated 

(7 = o.oir). 



pure two- level atoms [HI, [TH, Fig. [5] presents the 

noise spectra obtained for both quadratures in the case 
of resonant atomic excitation (A = 0) for a rather thin 
sample (C = 1), for different values of the reduced Rabi 
frequency U, r . Excess noise, with maximum value at 
zero frequency, is obtained for the phase quadrature for 
all field intensities. The amplitude quadrature presents 
squeezing centered at zero frequency for small light field 
amplitude (f2 r < 0.5F). As the Rabi frequency increases, 
the maximum squeezing shifts to increasing nonzero fre- 
quencies and disappears as f2 r > V. For large £l r , the 
intensity fluctuations present an excess noise peak cen- 
tered at the actual Rabi frequency f2o- Notice that, since 
fi r is the reduced Rabi frequency, the actual Rabi fre- 
quency associated to the specific two-level transition is 
^0 = (Qoo,ge)^r = T^f^r- As discussed in [lU, ampli- 
tude squeezing occurs in this case only in the limit of an 
optically thin medium, where the amount of squeezing 
is linearly dependent on C. However, at zero detuning, 
the mean field absorption, which is not taken into ac- 
count in our treatment, is considerable and also linearly 
dependent on C. Consequently, absorption of the inci- 
dent field should prevent the observation of squeezing 
in this regime. To avoid the effect of field absorption, 
one should consider situations in which the generalized 
Rabi frequency f2 = (A 2 + SI 2 ,) 1 / 2 is large. Such case 
corresponds to well resolved levels in the dressed-atom 
picture (38|. Noise spectra calculated for A = fl r = 1QT 
are presented in Fig. [3l for different values of the coop- 
erativity parameter. C = 1 corresponds to a thin atomic 
medium. The noise spectrum is essentially dominated 
by a peak occurring at w = f2. Squeezing occurs for 



the phase quadrature and, correspondingly, excess noise 
is present in the amplitude quadrature. The observed 
squeezing is due to four-wave mixing between the mean 
field at the carrier frequency ujl and the noise sidebands 
at frequencies u>l ± u> [l4| . In the dressed atom picture 
of the atoms driven by the incident mean field, a double 
A scheme occurs, involving the absorption (emission) of 
two driving field photons and the emission (absorption) 
of a photon from each of the two sidebands (see inset in 
Fig. [5]). This mechanism results in the buildup of corre- 
lations between the two sideband fluctuations, leading to 
quadrature squeezing . Such process is fully resonant 
with the energy levels of the dressed-atom when u> = O. 

As the medium's optical thickness increases, the noise 
feature around lu = broadens considerably and 
presents oscillations which are due to the frequency de- 
pendent phase mismatch between the carrier mean field 
and the noise sidebands [T3|. For large C broadband 
squeezing is present. More than 20% squeezing is ob- 
tained for C = 100. 

We notice from Fig. [3] that, except for the largest 
value of the cooperativity parameter C, there is no 
significant variation of the noise power around zero 
frequency. This behavior corresponds to closed two-level 
systems. Open two-level systems show an increase of 
light fluctuations at low frequency. 

The relative effect on the transmitted light fluctua- 
tions of the semiclassical atomic response and that of 
the atomic quantum fluctuations is illustrated in Fig. [4] 
The total spectrum for the phase quadrature noise and 
the contributions corresponding to the two terms on the 
right-hand side of Eq. [TH for A = 10r = fl r = 10T 
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FIG. 4: (Color online) Transition F g = -> F e = 1. Solid: to- 
tal noise spectrum for the phase quadrature (A = Q. r = 10F, 
7 = 0.01F, C — 100). Dashed: semiclassical contribution 
(fist term on the r.h.s. of Eq. [TSJ. Dotted: contribution of 
the atomic quantum fluctuations (second term on the r.h.s. 
of Eq. I18p . Inset: Dressed-atom picture for the four-wave- 
mixing process responsible for squeezing in two-level atoms. 
Hollow arrows: laser mean field; solid arrows: noise sidebands. 



and C — 100, are shown. As expected, the squeezing 
is only due to the semiclassical term contribution. For 
the chosen parameters, the atomic quantum fluctua- 
tions introduce a rather broad and smooth noise increase. 

The fact that, for this level scheme, squeezing occurs 
for one quadrature of the x polarized field, while the vac- 
uum fluctuations incident along the y polarization are un- 
affected, implies entanglement between two orthogonally 
polarized modes (39j . This can be explicitly verified by 
noticing that squeezing of the phase quadrature of field 1 
together with vacuum fluctuations in field 2 correspond 
to: 



a i 



A(a 2 



^ 2 



(31) 



Let now consider a + = ai +? 2 and 



— — 1 1 i j ■ i. ! f .._ — — ~ — J— ~ tllC 

field operators corresponding to the two orthogonal linear 
polarization at 45° with respect to the x and y axis. Then 

the operators X^ = and = a>± i J^ = +> — ) 

are two pairs of conjugate hermitian operators satisfying 
[X„,X V ] = 0, [Y^Yv] = 0, [X„,Y V ] = iS„ u . For these 
operators the inequality in Eq. 1311 becomes: 



A(Y+ + YJ) 2 + A(X + - A_) 2 < 2 



(32) 



which is sufficient to demonstrate continuous variable en- 
tanglement of the + and — polarization fields [231 ] - 



FIG. 5: (Color online) Transition F g = 1 -> F e = (the level 
scheme is shown using x as quantization axis; hollow arrow: 
driving field; solid arrows: spontaneous emission channels into 
the y polarized mode). Noise spectra of the amplitude (dot- 
ted) and phase (dashed) quadratures of the transmitted field 
with the same linear polarization as the driving field. The 
solid line corresponds to the noise of both quadratures of the 
output field with polarization orthogonal to the driving field 

(A = n r = lor, 7 = o.oir, c = 10). 



B. F a 



F e =0 



Excitation of a transition from a ground state with 
F g = 1 to a F e = excited state with a linearly (x) po- 
larized field corresponds to the coupling of the incident 
field to an open two-level system. Choosing the quan- 
tization axis along x determines that the incident field 
is coupled to the \F g = l,m g = 0) to \F e = 0,m e = 0) 
transition. The excited state can decay into either of 
the ground state Zeeman sublevels. Alternatively (in a 
different basis for the ground state manifold) , the config- 
uration corresponds to the excitation of one branch of a A 
system while vacuum field is acting on the second branch. 
A similar situation was studied in (33[, in the case of a 
resonant excitation (A = 0). Fig. [5]shows the noise spec- 
tra for both quadratures of the output fields with x and 
y polarizations. The spectra corresponding to the x po- 
larization are very similar to those obtained for a closed 
two-level system (Fig. [3]), except for the noise increase 
around zero frequency. Such a low-frequency feature in 
the noise spectra occurs in open transitions owing to the 
fluctuations introduced by spontaneous decay out of the 
two-level system. In addition, an increase in the noise 
power above the vacuum fluctuations is seen for the y 
polarized field. Both quadratures with this polarization 
experience the same noise contribution due to sponta- 
neous emission from the excited state. As a consequence 
of the light shift of the excited state produced by the driv- 
ing field, the noise for the y polarization peaks at lu =/= 0. 
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FIG. 6: (Color online) Transition F g = 1/2 F e = 1/2 
(the level scheme is shown using a; as quantization axis; hol- 
low arrow: driving field; solid arrows: spontaneous emission 
channels into the y polarized mode). Resonant noise spectra 
of the amplitude (solid) and phase (dashed) quadratures of 
the transmitted field with the same linear polarization as the 
driving field. Dotted: amplitude quadrature noise spectrum 
of the polarization perpendicular to the driving field, the spec- 
trum of the phase quadrature for this polarization coincides 
with the dashed line (A = 0, Q. r = lOf, 7 = O.Oir, C = 10). 



In addition to this low frequency peak, a much smaller 
local maximum occurs for w^!], almost invisible on the 
scale of Fig. O 



FIG. 7: (Color online) Transition F g = 1/2 -> F e = 1/2, 
A = 1GT. Amplitude (solid) and phase (dashed) quadratures 
of the transmitted field with the same linear polarization as 
the driving field. Amplitude (dash-dot) and phase (dotted) 
quadratures of the transmitted field with polarization perpen- 
dicular to the driving field. Inset: expanded low- frequency 
range (fi r = 10T, 7 = O.Oir, C = 10). 



tude (phase) noise and the y polarization phase (ampli- 
tude) noise. This approximate symmetry is broken for 
low noise frequencies (see inset on Fig. [7J|. Squeezing 
due to PSR occurs in this system for nonzero detuning. 
In addition, the inequality [21] is verified for u> ~ f2, giving 
rise to entanglement. 



C. F g = 1/2 F e = 1/2 

This configuration was analyzed in [20l | as a model 
system for the study of squeezing via PSR. Fig. [6] 
presents the corresponding spectra for resonant excita- 
tion (A = 0). The noise spectra for the driving field po- 
larization are analogous to those in Fig. FJ]for a two-level 
transition. In fact, for a choice of the quantization axis 
along x, the present situation corresponds to two two- 
level transitions coupled through spontaneous emission. 
Spontaneous emission is also responsible for the injec- 
tion of field fluctuations for the y polarization. Owing 
to the large symmetry of this configuration, the ampli- 
tude noise spectrum for the y polarization is equal to the 
phase noise spectrum of the x polarization. The phase 
noise spectrum for the y polarization presents enhanced 
noise if compared to the amplitude quadrature of the x 
polarization. 

The spectra for a non-resonant excitation (A = lOT) 
are presented in Fig. [7] For u> > T the spectra re- 
produce the features observed for non-resonant excita- 
tion of two- level systems (see Fig. [3]). Nearly identical 
noise spectra are obtained for the x polarization ampli- 



D. F g = 1 -> F e = 2 

We now consider a transition of the type F g > — > 
F e = F g + 1. For this class of transitions, no dark state 
exists within the ground level. At two-photon Raman 
resonance between ground state Zeeman sublevels, co- 
herence resonances occur that correspond to enhanced 
absorption, that is electromagnetically induced absorp- 
tion (EIA) [13, EH. The absorption spectra of a tran- 
sition of this typ e driven by a strong classical field were 
examined in [42j . If a driving field of circular polarization 
is used, owing to optical pumping, the system approxi- 
mates a pure two-level system. However, if linear polar- 
ization is used for the driving field, while the system is 
probed along the orthogonal polarization, the absorption 
spectrum presents a rather complex structure as a result 
of the different light shifts experienced by the Zeeman 
sublevels. A simple picture of this effect in terms of the 
dressed atom model is presented in [42j . 

The transmitted field noise spectra for the transition 
F g = 1 — > F e — 2 are presented in Fig. [5] for nonzero de- 
tuning (A = lOT) and large Rabi frequency (Sl r = 40T). 
A thick optical medium is considered (C = 100). For 
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FIG. 8: (Color online) Transition F g = 1 -> F e = 2 (the 
level scheme is shown using x as quantization axis; hollow ar- 
row: driving field; solid arrows: spontaneous emission chan- 
nels into the y polarized mode). Non- resonant noise spectra 
(A = 10r). Amplitude (solid) and phase (dashed) quadra- 
tures of the transmitted field with the same linear polarization 
as the driving field. Amplitude (dash-dot) and phase (dotted) 
quadratures of the transmitted field with polarization perpen- 
dicular to the driving field. Inset: expanded low-frequency 
range (Q r = 40I\ 7 = O.Oir, C = 100). 



lu > r, the spectra corresponding to the driving field po- 
larization are similar to the pure two-level spectra shown 
in Fig. [3J Unlike the pure two-level system, noise is in- 
troduced by spontaneous emission into the y polarized 
field. However, PSR squeezing is nevertheless present on 
the amplitude quadrature for this polarization. In ad- 
dition, the inequality in Eq. [3l] is verified resulting in 
entanglement. 

At low frequencies (uj < T), the noise spectra present 
features that are specific of such a multilevel system. Ex- 
cess noise peaking at uj — occurs for the two quadra- 
tures of the driving field polarization. The orthogo- 
nal polarization presents features which are peaked at 
uj > owing to the different Zeeman sublevel light shifts 
[42I |43| . Notice the significant squeezing near zero fre- 
quency of the phase quadrature. 



E. Effect of atomic motion 

We will now illustrate the effect of atomic motion on 
the noise spectra for an atomic sample with a Maxwell- 
Boltzmann velocity distribution. Having in mind the case 
of rubidium vapor at room temperature, we have used a 
Doppler width (FWHM) T Dopp = 90T. We consider as 
an example the F g = 1 — » F e = 2 transition (same as in 
Fig.©. 

Figure [5] shows the noise spectra for both field quadra- 



FIG. 9: (Color online) Transition F g = 1 -> F e = 2. Noise 
spectra for a Maxwell-Boltzmann atomic velocity distribu- 
tion and driving field tuned to resonance with zero velocity 
atoms. Amplitude (solid) and phase (dashed) quadratures 
of the transmitted field with the same linear polarization as 
the driving field. Amplitude (dash-dot) and phase (dotted) 
quadratures of the transmitted field with polarization perpen- 
dicular to the driving field. Inset: expanded low- frequency 
range (O r = 40r, 7 = O.Oir, C = 100, Y Dovp = 90r). 



tures in the case of a driving field resonant with the 
zero velocity atoms and a Rabi frequency il r — 40r for 
C = 100. Notice that, although the average detuning 
is zero in this case, the noise spectra are reminiscent of 
those obtained for atoms at rest with nonzero detuning. 
A local maximum (minimum)is observed for the ampli- 
tude (phase) quadrature noise of the driving field polar- 
ization for uj ~ f2o- The squeezing of the phase quadra- 
ture remains significant (~ 15%) in spite of the spreading 
of the field detuning caused by the Doppler effect. Notice 
that the squeezing of the amplitude quadrature noise of 
the y polarized field is not completely suppressed by the 
velocity distribution nor are the low frequency spectral 
features for both polarizations. 



IV. CONCLUSIONS 

We have developed a theoretical model supporting the 
numerical calculation of the noise spectra of light travers- 
ing an optically thick atomic medium, under a wide range 
of conditions. The model takes into account the com- 
plete Zeeman sublevel structure of the atomic transition 
and consequently the arbitrary polarization of the light 
field. In this paper, we have restricted our analysis to the 
case of a linearly polarized driving field while considering 
that the Zeeman sublevels are degenerate in the absence 
of light (zero magnetic field). We have calculated the 
noise spectra for the amplitude and phase quadratures of 
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the transmitted light with the same and the orthogonal 
polarization relative to the driving field. 

For atoms at rest, as can be produced in magneto- 
optical traps, the dominant feature of the transmitted 
field fluctuation spectra can be traced back to pure 
two- level atom effects [l3|. For nonzero optical detun- 
ing, excess noise is introduced in the amplitude quadra- 
ture while squeezing occurs for the phase quadrature. 
The noise variations are maximum for u> = Q, as a 
consequence of resonant four-wave mixing of the mean 
field and noise sidebands with the dressed-atom level 
structure. The propagation through the optically thick 
medium results in broadening and oscillatory structure 
of the noise spectrum owing to phase mismatch between 
the carrier and the fluctuation sidebands. Up to 30% 
squeezing occur for the parameter values used in the sim- 
ulations, chosen in order to correspond to typical experi- 
mental conditions, involving CW lasers and alkali atoms. 
The figure of C = 100 used in several of the present cal- 
culations corresponds to a realistic figure for magneto- 
optically trapped atom clouds. 

In addition to pure two-level atom features, novel ef- 
fects are present as a consequence of the Zeeman sub- 
level structure. Except when F g — 0, spontaneous emis- 
sion introduces noise into the polarization orthogonal to 
the driving field. Notwithstanding the random field fluc- 
tuations caused by spontaneous emission, nonlinear in- 
teraction between noise sidebands on both polarizations, 
may result (as in Figs. [7] and [8|) in vacuum squeezing 
for the orthogonal polarization. For the conditions used 
in our calculations, the noise introduced by the atomic 
quantum fluctuations does not completely mask the or- 
thogonal polarization vacuum squeezing. Furthermore, 
for transitions involving multiple Zeeman sublevels, new 
structures appear in the noise spectra at low frequency. 
They come from the different Stark shifts experienced by 
the Zeeman sublevels in the presence of the driving field. 
Such low frequency structures are reported here for the 
first time. 

The study presented here constitutes a necessary step 
towards a satisfactory understanding and control of 
atom-light interaction at the quantum level. This control 
can be applied to matter-light interfaces for quantum in- 
formation purposes, as the entanglement predicted here 
implies. Further developments may include the influence 



of nearby transitions as well as that of the spatial struc- 
ture of the light mode and atomic sample. Pulse propa- 
gation effects should also be considered. Such theoretical 
approach needs to be complemented with experimental 
tests, in particular using cold atomic samples. Work in 
this direction is currently underway. 



APPENDIX A: CALCULATION OF THE 
DIFFUSION COEFFICIENTS MATRIX 

The Heisenberg-Langevin equation (Eq. [5]) can be 
written in the form: 

^ = Z> (Ala) 

2?(<7y) EE 5^-4y,fci<7- H +7<4 . ( Alb ) 
kl 

where A = {Aij t ki} is the matrix appearing in Eq. 

[TUl The Langevin forces obey: ^/y (z,t) (z! ,t')\ = 

jj2D ljM 5 (z' - z") S(t- t'). The diffusion matrix D = 
{Dij : ki} can be calculated with the help of the general- 
ized Einstein relation [3?], Hfl : 

2Ai,« = (V (vi^l) - V {a l3 ) 4 - tryV (4) ) . 

(A2) 

From Eq. lAlb . by making use of crLcrju = \i) (J\ 
10 ( k \ = °\k 5 3l and Aj,mn = A*^ nm we get: 

= (A3) 
(v (o^) °i) = Aij,km (<rim) + Wij (o7fc) 

m 

(a i:j V (ah) ) = X] Alk < mi (°" m ^ + 7 ^ ^ ■ 

m 

The first of Eqs. IA3I results from the stationarity of the 
system. The mean values (<7ij) are obtained from the 
steady state solution of Eq. lAla . 
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